Four-Dimensional Mesoscale Liquid Model of Nucleus Resolves Chromatin's Radial Organization.

Recent advances chromatin capture, imaging techniques, and polymer modeling have dramatically enhanced quantitative understanding of chromosomal folding. However, the dynamism inherent in genome architectures due to physical and biochemical forces and their impact on nuclear architecture and cellular functions remains elusive. While imaging of chromatin in four dimensions is becoming more common, there is a conspicuous lack of physics-based computational tools appropriate for revealing the forces that shape nuclear architecture and dynamics. To this end, we have developed a multiphase liquid model of the nucleus, which can resolve chromosomal territories, compartments, and nuclear lamina using a physics-based and data-informed free-energy function. The model enables rapid hypothesis-driven prototyping of nuclear dynamics in four dimensions, thereby facilitating comparison with whole nucleus imaging experiments. As an application, we model the Drosophila nucleus and map phase diagram of various possible nuclear morphologies. We shed light on the interplay of adhesive and cohesive interactions which give rise to distinct radial organization seen in conventional, inverted, and senescent nuclear architectures. The results also show the highly dynamic nature of the radial organization, the disruption of which leads to significant variability in domain coarsening dynamics and consequently variability of chromatin architecture. The model also highlights the impact of oblate nuclear geometry and heterochromatin-subtype interactions on the global chromatin architecture and local asymmetry of chromatin compartments.

Chromatin alternates between A and B compartments at kilobase scale for subgenic organization.

Nuclear compartments are prominent features of 3D chromatin organization, but sequencing depth limitations have impeded investigation at ultra fine-scale. CTCF loops are generally studied at a finer scale, but the impact of looping on proximal interactions remains enigmatic. Here, we critically examine nuclear compartments and CTCF loop-proximal interactions using a combination of in situ Hi-C at unparalleled depth, algorithm development, and biophysical modeling. Producing a large Hi-C map with 33 billion contacts in conjunction with an algorithm for performing principal component analysis on sparse, super massive matrices (POSSUMM), we resolve compartments to 500 bp. Our results demonstrate that essentially all active promoters and distal enhancers localize in the A compartment, even when flanking sequences do not. Furthermore, we find that the TSS and TTS of paused genes are often segregated into separate compartments. We then identify diffuse interactions that radiate from CTCF loop anchors, which correlate with strong enhancer-promoter interactions and proximal transcription. We also find that these diffuse interactions depend on CTCF's RNA binding domains. In this work, we demonstrate features of fine-scale chromatin organization consistent with a revised model in which compartments are more precise than commonly thought while CTCF loops are more protracted.

Structural reorganization and relaxation dynamics of axially stressed chromosomes.

Chromosomes endure mechanical stresses throughout the cell cycle; for example, resulting from the pulling of chromosomes by spindle fibers during mitosis or deformation of the nucleus during cell migration. The response to physical stress is closely related to chromosome structure and function. Micromechanical studies of mitotic chromosomes have revealed them to be remarkably extensible objects and informed early models of mitotic chromosome organization. We use a data-driven, coarse-grained polymer modeling approach to explore the relationship between the spatial organization of individual chromosomes and their emergent mechanical properties. In particular, we investigate the mechanical properties of our model chromosomes by axially stretching them. Simulated stretching led to a linear force-extension curve for small strain, with mitotic chromosomes behaving about 10-fold stiffer than interphase chromosomes. Studying their relaxation dynamics, we found that chromosomes are viscoelastic solids with a highly liquid-like, viscous behavior in interphase that becomes solid-like in mitosis. This emergent mechanical stiffness originates from lengthwise compaction, an effective potential capturing the activity of loop-extruding SMC complexes. Chromosomes denature under large strains via unraveling, which is characterized by opening of large-scale folding patterns. By quantifying the effect of mechanical perturbations on the chromosome's structural features, our model provides a nuanced understanding of in vivo mechanics of chromosomes.

An associative memory Hamiltonian model for DNA and nucleosomes.

A model for DNA and nucleosomes is introduced with the goal of studying chromosomes from a single base level all the way to higher-order chromatin structures. This model, dubbed the Widely Editable Chromatin Model (WEChroM), reproduces the complex mechanics of the double helix including its bending persistence length and twisting persistence length, and the temperature dependence of the former. The WEChroM Hamiltonian is composed of chain connectivity, steric interactions, and associative memory terms representing all remaining interactions leading to the structure, dynamics, and mechanical characteristics of the B-DNA. Several applications of this model are discussed to demonstrate its applicability. WEChroM is used to investigate the behavior of circular DNA in the presence of positive and negative supercoiling. We show that it recapitulates the formation of plectonemes and of structural defects that relax mechanical stress. The model spontaneously manifests an asymmetric behavior with respect to positive or negative supercoiling, similar to what was previously observed in experiments. Additionally, we show that the associative memory Hamiltonian is also capable of reproducing the free energy of partial DNA unwrapping from nucleosomes. WEChroM is designed to emulate the continuously variable mechanical properties of the 10nm fiber and, by virtue of its simplicity, is ready to be scaled up to molecular systems large enough to investigate the structural ensembles of genes. WEChroM is implemented in the OpenMM simulation toolkits and is freely available for public use.

Interphase chromosomes of the Aedes aegypti mosquito are liquid crystalline and can sense mechanical cues.

We use data-driven physical simulations to study the three-dimensional architecture of the Aedes aegypti genome. Hi-C maps exhibit both a broad diagonal and compartmentalization with telomeres and centromeres clustering together. Physical modeling reveals that these observations correspond to an ensemble of 3D chromosomal structures that are folded over and partially condensed. Clustering of the centromeres and telomeres near the nuclear lamina appears to be a necessary condition for the formation of the observed structures. Further analysis of the mechanical properties of the genome reveals that the chromosomes of Aedes aegypti, by virtue of their atypical structural organization, are highly sensitive to the deformation of the nuclei. This last finding provides a possible physical mechanism linking mechanical cues to gene regulation.

Shaping the genome via lengthwise compaction, phase separation, and lamina adhesion.

The link between genomic structure and biological function is yet to be consolidated, it is, however, clear that physical manipulation of the genome, driven by the activity of a variety of proteins, is a crucial step. To understand the consequences of the physical forces underlying genome organization, we build a coarse-grained polymer model of the genome, featuring three fundamentally distinct classes of interactions: lengthwise compaction, i.e., compaction of chromosomes along its contour, self-adhesion among epigenetically similar genomic segments, and adhesion of chromosome segments to the nuclear envelope or lamina. We postulate that these three types of interactions sufficiently represent the concerted action of the different proteins organizing the genome architecture and show that an interplay among these interactions can recapitulate the architectural variants observed across the tree of life. The model elucidates how an interplay of forces arising from the three classes of genomic interactions can drive drastic, yet predictable, changes in the global genome architecture, and makes testable predictions. We posit that precise control over these interactions in vivo is key to the regulation of genome architecture.

The interplay of chromatin phase separation and lamina interactions in nuclear organization.

The genetic material of eukaryotes is segregated into transcriptionally active euchromatin and silent heterochromatin compartments. The spatial arrangement of chromatin compartments evolves over the course of cellular life in a process that remains poorly understood. The latest nuclear imaging experiments reveal a number of dynamical signatures of chromatin that are reminiscent of active multiphase liquids. This includes the observations of viscoelastic response, coherent motions, Ostwald ripening, and coalescence of chromatin compartments. There is also growing evidence that liquid-liquid phase separation of protein and nucleic acid components is the underlying mechanism for the dynamical behavior of chromatin. To dissect the organizational and dynamical implications of chromatin's liquid behavior, we have devised a phenomenological field-theoretic model of the nucleus as a multiphase condensate of liquid chromatin types. Employing the liquid chromatin model of the Drosophila nucleus, we have carried out an extensive set of simulations with an objective to shed light on the dynamics and chromatin patterning observed in the latest nuclear imaging experiments. Our simulations reveal the emergence of experimentally detected mesoscale chromatin channels and spheroidal droplets which arise from the dynamic interplay of chromatin type to type interactions and intermingling of chromosomal territories. We also quantitatively reproduce coherent motions of chromatin domains observed in displacement correlation spectroscopy measurements which are explained within the framework of our model by phase separation of chromatin types operating within constrained intrachromosomal and interchromosomal boundaries. Finally, we illuminate the role of heterochromatin-lamina interactions in the nuclear organization by showing that these interactions enhance the mobility of euchromatin and indirectly introduce correlated motions of heterochromatin droplets.

A Liquid State Perspective on Dynamics of Chromatin Compartments.

The interior of the eukaryotic cell nucleus has a crowded and heterogeneous environment packed with chromatin polymers, regulatory proteins, and RNA molecules. Chromatin polymer, assisted by epigenetic modifications, protein and RNA binders, forms multi-scale compartments which help regulate genes in response to cellular signals. Furthermore, chromatin compartments are dynamic and tend to evolve in size and composition in ways that are not fully understood. The latest super-resolution imaging experiments have revealed a much more dynamic and stochastic nature of chromatin compartments than was appreciated before. An emerging mechanism explaining chromatin compartmentalization dynamics is the phase separation of protein and nucleic acids into membraneless liquid condensates. Consequently, concepts and ideas from soft matter and polymer systems have been rapidly entering the lexicon of cell biology. In this respect, the role of computational models is crucial for establishing a rigorous and quantitative foundation for the new concepts and disentangling the complex interplay of forces that contribute to the emergent patterns of chromatin dynamics and organization. Several multi-scale models have emerged to address various aspects of chromatin dynamics, ranging from equilibrium polymer simulations, hybrid non-equilibrium simulations coupling protein binding and chromatin folding, and mesoscopic field-theoretic models. Here, we review these emerging theoretical paradigms and computational models with a particular focus on chromatin's phase separation and liquid-like properties as a basis for nuclear organization and dynamics.

The Nucleome Data Bank: web-based resources to simulate and analyze the three-dimensional genome.

We introduce the Nucleome Data Bank (NDB), a web-based platform to simulate and analyze the three-dimensional (3D) organization of genomes. The NDB enables physics-based simulation of chromosomal structural dynamics through the MEGABASE + MiChroM computational pipeline. The input of the pipeline consists of epigenetic information sourced from the Encode database; the output consists of the trajectories of chromosomal motions that accurately predict Hi-C and fluorescence insitu hybridization data, as well as multiple observations of chromosomal dynamics in vivo. As an intermediate step, users can also generate chromosomal sub-compartment annotations directly from the same epigenetic input, without the use of any DNA-DNA proximity ligation data. Additionally, the NDB freely hosts both experimental and computational structural genomics data. Besides being able to perform their own genome simulations and download the hosted data, users can also analyze and visualize the same data through custom-designed web-based tools. In particular, the one-dimensional genetic and epigenetic data can be overlaid onto accurate 3D structures of chromosomes, to study the spatial distribution of genetic and epigenetic features. The NDB aims to be a shared resource to biologists, biophysicists and all genome scientists. The NDB is available at https://ndb.rice.edu.

Exploring chromosomal structural heterogeneity across multiple cell lines.

Using computer simulations, we generate cell-specific 3D chromosomal structures and compare them to recently published chromatin structures obtained through microscopy. We demonstrate using machine learning and polymer physics simulations that epigenetic information can be used to predict the structural ensembles of multiple human cell lines. Theory predicts that chromosome structures are fluid and can only be described by an ensemble, which is consistent with the observation that chromosomes exhibit no unique fold. Nevertheless, our analysis of both structures from simulation and microscopy reveals that short segments of chromatin make two-state transitions between closed conformations and open dumbbell conformations. Finally, we study the conformational changes associated with the switching of genomic compartments observed in human cell lines. The formation of genomic compartments resembles hydrophobic collapse in protein folding, with the aggregation of denser and predominantly inactive chromatin driving the positioning of active chromatin toward the surface of individual chromosomal territories.

Mesoscale Liquid Model of Chromatin Recapitulates Nuclear Order of Eukaryotes.

The nuclear envelope segregates the genome of Eukaryota from the cytoplasm. Within the nucleus, chromatin is further compartmentalized into architectures that change throughout the lifetime of the cell. Epigenetic patterns along the chromatin polymer strongly correlate with chromatin compartmentalization and, accordingly, also change during the cell life cycle and at differentiation. Recently, it has been suggested that subnuclear chromatin compartmentalization might result from a process of liquid-liquid phase separation orchestrated by the epigenetic marking and operated by proteins that bind to chromatin. Here, we translate these observations into a diffuse interface model of chromatin, which we named the mesoscale liquid model of nucleus. Using this streamlined continuum model of the genome, we study the large-scale rearrangements of chromatin that happen at different stages of the growth and senescence of the cell and during nuclear inversion events. In particular, we investigate the role of droplet diffusion, fluctuations, and heterochromatin-lamina interactions during nuclear remodeling. Our results indicate that the physical process of liquid-liquid phase separation, together with surface effects, is sufficient to recapitulate much of the large-scale morphology and dynamics of chromatin along the life cycle of cells.

Walking along chromosomes with super-resolution imaging, contact maps, and integrative modeling.

Chromosome organization is crucial for genome function. Here, we present a method for visualizing chromosomal DNA at super-resolution and then integrating Hi-C data to produce three-dimensional models of chromosome organization. Using the super-resolution microscopy methods of OligoSTORM and OligoDNA-PAINT, we trace 8 megabases of human chromosome 19, visualizing structures ranging in size from a few kilobases to over a megabase. Focusing on chromosomal regions that contribute to compartments, we discover distinct structures that, in spite of considerable variability, can predict whether such regions correspond to active (A-type) or inactive (B-type) compartments. Imaging through the depths of entire nuclei, we capture pairs of homologous regions in diploid cells, obtaining evidence that maternal and paternal homologous regions can be differentially organized. Finally, using restraint-based modeling to integrate imaging and Hi-C data, we implement a method-integrative modeling of genomic regions (IMGR)-to increase the genomic resolution of our traces to 10 kb.

Deciphering the structure of the condensin protein complex.

Protein assemblies consisting of structural maintenance of chromosomes (SMC) and kleisin subunits are essential for the process of chromosome segregation across all domains of life. Prokaryotic condensin belonging to this class of protein complexes is composed of a homodimer of SMC that associates with a kleisin protein subunit called ScpA. While limited structural data exist for the proteins that comprise the (SMC)-kleisin complex, the complete structure of the entire complex remains unknown. Using an integrative approach combining both crystallographic data and coevolutionary information, we predict an atomic-scale structure of the whole condensin complex, which our results indicate being composed of a single ring. Coupling coevolutionary information with molecular-dynamics simulations, we study the interaction surfaces between the subunits and examine the plausibility of alternative stoichiometries of the complex. Our analysis also reveals several additional configurational states of the condensin hinge domain and the SMC-kleisin interaction domains, which are likely involved with the functional opening and closing of the condensin ring. This study provides the foundation for future investigations of the structure-function relationship of the various SMC-kleisin protein complexes at atomic resolution.

Anomalous diffusion, spatial coherence, and viscoelasticity from the energy landscape of human chromosomes.

The nucleus of a eukaryotic cell is a nonequilibrium system where chromatin is subjected to active processes that continuously rearrange it over the cell's life cycle. Tracking the motion of chromosomal loci provides information about the organization of the genome and the physical processes shaping that organization. Optical experiments report that loci move with subdiffusive dynamics and that there is spatially coherent motion of the chromatin. We recently showed that it is possible to predict the 3D architecture of genomes through a physical model for chromosomes that accounts for the biochemical interactions mediated by proteins and regulated by epigenetic markers through a transferable energy landscape. Here, we study the temporal dynamics generated by this quasi-equilibrium energy landscape assuming Langevin dynamics at an effective temperature. Using molecular dynamics simulations of two interacting human chromosomes, we show that the very same interactions that account for genome architecture naturally reproduce the spatial coherence, viscoelasticity, and the subdiffusive behavior of the motion in interphase chromosomes as observed in numerous experiments. The agreement between theory and experiments suggests that even if active processes are involved, an effective quasi-equilibrium landscape model can largely mimic their dynamical effects.

The Energetics and Physiological Impact of Cohesin Extrusion.

Cohesin extrusion is thought to play a central role in establishing the architecture of mammalian genomes. However, extrusion has not been visualized in vivo, and thus, its functional impact and energetics are unknown. Using ultra-deep Hi-C, we show that loop domains form by a process that requires cohesin ATPases. Once formed, however, loops and compartments are maintained for hours without energy input. Strikingly, without ATP, we observe the emergence of hundreds of CTCF-independent loops that link regulatory DNA. We also identify architectural "stripes," where a loop anchor interacts with entire domains at high frequency. Stripes often tether super-enhancers to cognate promoters, and in B cells, they facilitate Igh transcription and recombination. Stripe anchors represent major hotspots for topoisomerase-mediated lesions, which promote chromosomal translocations and cancer. In plasmacytomas, stripes can deregulate Igh-translocated oncogenes. We propose that higher organisms have coopted cohesin extrusion to enhance transcription and recombination, with implications for tumor development.

De novo prediction of human chromosome structures: Epigenetic marking patterns encode genome architecture.

Inside the cell nucleus, genomes fold into organized structures that are characteristic of cell type. Here, we show that this chromatin architecture can be predicted de novo using epigenetic data derived from chromatin immunoprecipitation-sequencing (ChIP-Seq). We exploit the idea that chromosomes encode a 1D sequence of chromatin structural types. Interactions between these chromatin types determine the 3D structural ensemble of chromosomes through a process similar to phase separation. First, a neural network is used to infer the relation between the epigenetic marks present at a locus, as assayed by ChIP-Seq, and the genomic compartment in which those loci reside, as measured by DNA-DNA proximity ligation (Hi-C). Next, types inferred from this neural network are used as an input to an energy landscape model for chromatin organization [Minimal Chromatin Model (MiChroM)] to generate an ensemble of 3D chromosome conformations at a resolution of 50 kilobases (kb). After training the model, dubbed Maximum Entropy Genomic Annotation from Biomarkers Associated to Structural Ensembles (MEGABASE), on odd-numbered chromosomes, we predict the sequences of chromatin types and the subsequent 3D conformational ensembles for the even chromosomes. We validate these structural ensembles by using ChIP-Seq tracks alone to predict Hi-C maps, as well as distances measured using 3D fluorescence in situ hybridization (FISH) experiments. Both sets of experiments support the hypothesis of phase separation being the driving process behind compartmentalization. These findings strongly suggest that epigenetic marking patterns encode sufficient information to determine the global architecture of chromosomes and that de novo structure prediction for whole genomes may be increasingly possible.

Transferable model for chromosome architecture.

In vivo, the human genome folds into a characteristic ensemble of 3D structures. The mechanism driving the folding process remains unknown. We report a theoretical model for chromatin (Minimal Chromatin Model) that explains the folding of interphase chromosomes and generates chromosome conformations consistent with experimental data. The energy landscape of the model was derived by using the maximum entropy principle and relies on two experimentally derived inputs: a classification of loci into chromatin types and a catalog of the positions of chromatin loops. First, we trained our energy function using the Hi-C contact map of chromosome 10 from human GM12878 lymphoblastoid cells. Then, we used the model to perform molecular dynamics simulations producing an ensemble of 3D structures for all GM12878 autosomes. Finally, we used these 3D structures to generate contact maps. We found that simulated contact maps closely agree with experimental results for all GM12878 autosomes. The ensemble of structures resulting from these simulations exhibited unknotted chromosomes, phase separation of chromatin types, and a tendency for open chromatin to lie at the periphery of chromosome territories.

A Stochastic Algorithm for the Isobaric-Isothermal Ensemble with Ewald Summations for All Long Range Forces.

We present an algorithm termed COMPEL (COnstant Molecular Pressure with Ewald sum for Long range forces) to conduct simulations in the NPT ensemble. The algorithm combines novel features recently proposed in the literature to obtain a highly efficient and accurate numerical integrator. COMPEL exploits the concepts of molecular pressure, rapid stochastic relaxation to equilibrium, exact calculation of the contribution to the pressure of long-range nonbonded forces with Ewald summation, and the use of Trotter expansion to generate a robust, highly stable, symmetric, and accurate algorithm. Explicit implementation in the MOIL program and illustrative numerical examples are discussed.